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ABSTRACT 

o\ 

ON \ The genus statistics of isodensity contours has become a well-established tool in cosmology. In 

this Letter we place the genus in the wider framework of a complete family of morphological 
descriptors. These are known as the Minkowski functionals, and we here apply them for the first 
time to isodensity contours of a continuous random field. By taking two equivalent approaches, 
one through differential geometry, the other through integral geometry, we derive two comple- 
mentary formulae suitable for numerically calculating the Minkowski functionals. As an example 
we apply them to simulated Gaussian random fields and compare the outcome to the analytically 
known results, demonstrating that both are indeed well suited for numerical evaluation. The 
code used for calculating all Minkowski functionals is available from the authors. 

m 

£N| ' Subject headings: large-scale structure of universe, methods: statistical 

O" 

~~^5 1. Introduction 

The ambitious goal of characterizing the morphology of cosmic structure in a unique and complete way 
leads to the field of integral geometry. There, this problem has been well-posed and Hadwiger's theorem 



( Hadwiger 1957 , see also Klain 1995 for a proof) gives a powerful answer: In d spatial dimensions the global 



morphological properties (defined as those which satisfy motional invariance and additivity) of any pattern 

•i-H 

X 



can be completely characterized by d + 1 numbers, the so-called Minkowski functionals (Minkowski 1903). 
Furthermore, in the interesting cases of two and three dimensions, all the Minkowski functionals can be 
jJJ ' intuitively interpreted in terms of well-known geometric quantities. 

& . 

One of them is the Euler characteristic \i or equivalently the genus g — 1 — \, which has long been used as 

a statistical tool in cosmology. The standard analysis technique starts from a density or temperature field. 

For a point process it is usually constructed by applying a smoothing filter with a constant width larger 

than the mean separation of points. For this field, the Euler characteristic of an isocontour at level v is 

obtained, according to the Gauss-Bonnet theorem, through a surface integration of the Gaussian curvature. 

This threshold level can be employed as a diagnostic parameter, and the genus represented as a function 



of v (see Weinberg et al. 1987 , Mclott 199C , Coles et al. 1996 and references therein) can be compared 



with analytically known values for Gaussian random fields ( Doroshkevich 1970 ) and even for some nonlinear 
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dynamical models such as solutions of the Lagrangian perturbation theory (Matsubara & Suto 1996; Seto 



et al. 1997) 



This method for calculating the genus or, as advocated in the present Letter, all four Minkowski functionals, 



suffers from erasing small-scale information (ter Haar Romeny et al. 1991) when going from point sets 
to smooth density fields. While this may actually be an advantage when trying to deemphasize poorly 
understood effects such as nonlinear evolution and biasing, some discriminatory power is lost. 



In this situation, Mecke et al. (1994) introduced Minkowski functionals of point processes into cosmology. 
Proceeding as directly as possible they decorate each point with a ball of radius R, which becomes the 
diagnostic parameter for this method of analysis. Thereby a union set of balls is constructed and probes 
the morphological properties of the point set. These are measured in a complete and unique way using 
Minkowski functionals, so powerful methods of integral geometry ( Weil 1983 ) become applicable. Most 
notably, Minkowski functionals contain a variety of other descriptors such as the "void probability function" 



(White 1979), and depend on the complete hierarchy of correlations (Stratonovich 1963; Mecke et al. 1994) 



Analytically known results for the Poisson process (Mecke & Wagner 1991) provide a useful standard of 
reference; the additivity property in integral geometry leads to robustness against various sources of error 



(Kerscher et al. 1997a); the sensitivity to small-scale clustering properties allows for efficient discrimination, 



as has been demonstrated in recent applications to the Abell/ACO cluster catalogue ( Kerscher et al. 19971 ) 



and the IRAS 1.2Jy catalogue ( Kerscher et al. 1997a ). A brief tutorial on Minkowski functionals in cosmology 



can be found in Schmalzing et al. (1995) 



In this Letter we advocate the use of the whole family of Minkowski functionals to analyze the continuous 
random fields encountered in cosmology. Thereby we obtain more information than by measuring the genus 



alone. Furthermore, recent work ( Kerscher et al. 1997a , see also Buchert 1995 for a direct comparison) has 
shown that other Minkowski functionals such as the surface area and the integrated mean curvature can 
possess greater discriminatory power than the Euler characteristic. Apart from this practical advantage, 



which was already pointed out for the area by Ryden et al. (1989), the use of all the Minkowski functionals 
leads to a conceptual generalization of genus statistics, incorporates it into the firm mathematical framework 
of integral geometry and stereology and allows for a straightforward evaluation, including exact correction 
for boundary effects. 

The Letter is organized as follows. Section 2 uses the methods of differential geometry to reduce the eval- 
uation of Minkowski functionals to a spatial average of Koenderink invariants. Section 3 employs integral 
geometry, namely Crofton's intersection formula, to derive a second calculational method that relies on the 
representation of the density field by its values at lattice points. Section 4 gives the analytical formulae for all 
the Minkowski functionals of Gaussian random fields, while Section 5 compares them to the results obtained 
by applying our two approximation formulae to numerically simulated Gaussian random fields. Section 6 
provides an outlook. 



Koenderink Invariants 



Consider a random field v(x) on a d-dimensional support V C IR . For a given threshold v the excursion 
set F v is the set of all points x with z/(x) > v. We wish to calculate its Minkowski functionals v k {v) per 
unit volume. The volume functional Vq is simply given by volume integration of a Heaviside step function 



0, normalized to the whole volume \T>\ 
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In three dimensions, the other Minkowski functionals correspond to quantities known from differential geom- 



etry, namely the surface area, the integrated mean curvature and the Euler characteristic (Schneider 1993) 
They can be evaluated by a surface integration 



Vfc(z/) 
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The local Minkowski functionals 
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are nothing but appropriately normalized invariants made from the inverse radii of curvature R\ and i?2 of 
the surface orientated towards lower density values. 



Using methods developed by Koenderink ( Kocnderink 1984 ; tcr Haar Romcny et al. 1991 ) in two-dimensional 
image analysis, it is possible to express these local curvatures in terms of geometric invariants formed from 
first and second derivatives ( Schmalzing 199q ). The surface integration (pj) is equivalent to taking the spatial 
mean over the whole volume, with a gradient for the surface element and a delta function for selecting the 
isodensity contour added as factors. The resulting formula is 



Vfc(^) 
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If the density field is sampled at some lattice points, we can therefore estimate the Minkowski functionals by 
calculating the derivatives at each lattice point and replacing first the (S-function with a bin of finite width 
and the average over all space with a sum over lattice points only. Analogous formulae can be derived in two 
dimensions, so the most promising cosmological application of this particular method is an elegant way of 
calculating Minkowski functionals for fields on the sphere, like all-sky maps of the CMB (for existing work, 
see Bmoot ct al. 1994 , Colley et al. 1996] and especially Torres ct al. 1995| ). 



Crofton's Intersection Formula 



Another powerful approach to morphology of density fields comes through integral geometry (Hadwiger 
1957 ; [Bantalo 1976 ). As above we consider a random field i/(x) and measure the Minkowski functionals 



77) 



(v) of its excursion set over a threshold v. 



To do this we use Crofton's formula ( Crofton 1868 ). For a body K in d dimensions we consider an arbitrary 
fc-dimensional hyperplane E and calculate the Euler characteristic x of the intersection KC\E in k dimensions. 
Integrating this quantity over the space Sj, of all conceivable hyperplanesf] we obtain the feth Minkowski 



5 The integration measure dfj, k (E) is normalized to give f (d) d[i k (E) = 1. 



functional of the body K in d dimensions. Crofton's formula states that[] 

v[ d) (K) = -^— /" d» k (E) X (k HKnE). (5) 



If the body K is the excursion set of a homogeneous and isotropic random field sampled at L points of a 
cubic lattice of spacing a, we can as well average over the set L k of all dual lattice hyperplanes only. The 
integral is then replaced by a summation over these hyperplanes with an additional normalization factor and 
we obtain 

k v" ' / . / r m 



vr>{K) = -?t- E 4r kJ ^ l ^(KnE). ( 6 ) 
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However, the Euler characteristic of a body in k dimensions can be approximated by counting the numbers 



Nj(K) of all elementary lattice cells of dimension j within the body (Adler 1981). In three dimensions, for 
example, N 3 (K) gives the number of cubes within the body, N 2 (K) and Ni(K) arc the numbers of faces 
and edges of the lattice cubes, respectively, while Nq(K) gives the number of lattice points contained in K. 
In general 

k 

x (k) (KnE) = ^(-iyNj(KnE). (7) 

3=0 

When summing equation (]7|) over all lattice hyperplanes, we notice that j-dimensional lattice cells of the 
bodies KP\E result from intersections of a (d — k + j) -dimensional lattice cell belonging to the body K 

k 

J 



with A:-dimensional hyperplanes E, so there are ( • -' I of them. So we can actually perform the 



summation 

E 



N t (KnE) =( K d ~ k J +J ^j N d - k+j (K) (8) 

and end up with the compact formula 



For a homogeneous and isotropic random field in three dimensions the Minkowski functionals are thus given 

by 

vi (10 = -^(-lN 3 (K) + IN 2 (K)) 

MK) = jyd^W - IN 2 (K) + lm(K)) 

MK) = -L(-N 3 (K) + N 2 {K) - N^K) + N (K)). 

CI _/v 

Note that the last relation for the Euler characteristic V3 = x/l^l P er un it volume is not obtained from 
Crofton's formula but is already given as a special case of equation (jTJ) which results from the Morse theorem 
(Morse fc Cairns 1969]). Being conceptually simpler than the other approximations, it was proved rigorously 



(10) 



by Adler (1981) and used by Coles et al. (1996) to calculate the genus of isodensity contours in three 



6 uJk is the volume of the unit ball in k dimensions, so luq = 1, ui\ = 2, u>2 = 1C and u>s = ^7 



dimensions. Adding Crofton's formula gives all the Minkowski functionals without further computational 
effort, since the calculation requires nothing but the counting of all elementary lattice cells, which is already 
necessary for the genus. 



Tomita's Formulae for Statistics of Random Interfaces 



Like the genus ( Doroshkcvich 1970| ), all Minkowski functionals are analytically known for a Gaussian random 
field i/(x) in d dimensions, due to a highly instructive article by Tomita (1990) . Again, consider the excursion 



set over a threshold v and calculate its mean Minkowski functionals vi {v) per unit volume. 

The result depends on only two parameters a and A, which in turn depend on the value of the correlation 
function £(0) and its second derivative £"(0) at zero through 



<t = £(0), A = ^"(0)1/(27^(0)). (11) 

Both are easily measured from a given realization of the random field, since £(0) = {v 2 ) is the variance of 



the field itself and |£"(0)| = (v t ) is the variance of any of its first derivatives (Adler 1981) 



With these two parameters all Minkowski functionals in arbitrary dimension d can be calculated for a 
Gaussian random field ( Tomita 1990 ). For the volume functional we obtain]] 
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while all other Minkowski functionals are given by 
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Putting d = 3 we obtain the Minkowski functionals in three dimensions^] 

vo(v) = 
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v 2 (^) = 
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Gaussian Random Fields 



Gaussian random fields provide a useful standard of reference and are always good for a pedagogical example. 
Here we use a scale-invariant power spectrum to construct Gaussian random fields on 128 3 points in a unit 



7 The function <3>(x) = -7= / dtexp(— i 2 ) denotes the Gaussian error integral, while H n (x) = (— -r-) I— 
defines a Hermite function of order n. 

8 To relieve the notation of the overall normalization y/a, we use the dimensionless argument u, with v = Uy/a. 



cxp(-ix 2 ) 



cube. We apply a Gaussian smoothing filter of width 16/128 = 0.125, which is chosen large in order to give 
a large variation between realizations. The field is normalized to (is) = and (y 2 ) = 1. Then, we calculate 
the numerical estimates and the analytical predictions at 100 threshold values between —4 and 4 and obtain 
the mean and variance over 200 such curves. 

Figure [0 shows three areas in each plot. One is measured by averaging Koenderink invariants as in equation 
(Q), the second by using Crofton's formula (|l0|). However, these two are completely indistinguishable, so the 
two approximations are indeed equivalent. 

A third area (shaded, dotted lines) is obtained by estimating the parameters a and A from each random 
field and computing the analytical mean Minkowski functionals according to (JIJ). Obviously the mea- 
sured mean values fit the analytical expectations extremely well; deviations are far from significant. Since 
the analytical curves incorporate only statistical fluctuations in the determination of the parameters, their 
variance is reduced compared to the measured curves, which contain the full variance between and within 
samples. Obviously, surface area vi and integrated mean curvature V2 vary almost exclusively through pa- 
rameter fluctuations, so these measures enhance the discriminatory power compared to considering the Eulcr 
characteristic V3 alone. 



6. Summary and Outlook 

We have extended the scope of genus statistics by considering all four Minkowski functionals in three dimen- 
sions. Two independent approaches gave us two different approximation formulae which we have tested by 
application to a Gaussian random field. Furthermore, we verified that our measurements fit the analytically 
known Minkowski functionals of a Gaussian random field. 

Even without analytical reference values, we have two complementary numerical formulae at our disposal 
which both rely on the approximation of a smooth density field by values at lattice points. The agreement 
between the two can be employed to determine the necessary smoothing length. Too little smoothing would 
be easily detectable by discrepancies between the two formulae. 

Forthcoming applications of the presented method will give intuitive interpretations of our novel access 



to cosmic morphology. Practical issues such as boundary corrections (sec Kcrschcr ct al. 1996 for point 
processes) , and further tests of the robustness and discriminatory power are necessary beyond the promising 
results given here. 
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Figure Captions 



Fig. 1. — The four Minkowski functionals of isodensity contours of a Gaussian random field in three 
dimensions. The central lines show the exact mean values while the area between top and bottom lines 
indicates the statistical variance of our results over 200 realizations. Actually there are three areas in each 
plot, corresponding to the three formulae (0), (nG) and (14). However, the two independent approximation 



formulae (no shading) perform so well that their results are identical. Also, the agreement with analytically 
predicted mean values (shaded) is remarkable. 
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